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The nature of the molecular gas system in the core of NGC 1275 
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ABSTRACT 

We present near-infrared integral field spectroscopy of the central kiloparsec of NGC 1275 at 
the heart of the Perseus cluster of galaxies, obtained with the UIST IFU on UKIRT. The nu- 
clear ro-vibrational H2 emission is spatially resolved and is likely to originate approximately 
50 pc from the active nucleus. The Paa emission is, by contrast, spatially unresolved. The 
requirements for thermal excitation of the H2 by nuclear X-radiation, its kinematics on sub- 
arcsec scales, and its stability against self-gravity, together suggest that the observed H2 is 
part of clumpy disk rotating about the radio-jet axis. The sharp jump in the H2 velocity across 
the nucleus implies a black hole mass of 3.4 x 10 s M Q , with a systematic error of ±0.18 dex 
due to the uncertainty in the radio jet inclination. This agrees well with the value implied by 
the empirical correlation between black hole mass and stellar velocity dispersion for nearby 
elliptical galaxies, and is ~ 100 times the stellar mass in this region. 

Key words: galaxies:active - galaxies: cooling flow - galaxiesdndividual: NGC 1275 - in- 
frared:galaxies 



1 INTRODUCTION 



For several years, NGC 1275 at the centre of the Perseus cluster 
held the unique distinction of being the only galaxy in a cooling 
flow cluster for which CO emission from molecular gas had been 
detected (Inoue et al. 1996; Bridges & Irwin 1998 and references 
therein). The fact that it is also a well-known galaxy merger (as 
well as having been classified a Seyfert 1, a BL Lac and an FR I 
radio source) served to obscure any causal connection between the 
CO and the cooling flow, so the failure to detect molecular gas in 
other such systems led to a heated debate over the validity of the 
cooling flow model. The latter model proposed that hot gas at the 
centres of relaxed clusters of galaxies was cooling out from the hot 
phase at rates of hundreds or thousands of solar masses per year 
(Fabian 1994). Recent developments on two separate fronts have, 
however, helped to bridge the impasse. 

Firstly, X-ray data from Chandra and XMM-Newton show that 
cooling rates were vastly overestimated in the past (e.g. Schmidt, 
Allen & Fabian 2001), and X-ray grating spectra also reveal a 
deficit of line emission from gas cooling below temperatures ~ 
2vMai/3 (Peterson et al. 2003). Explanations for these findings in- 
clude rapid mixing of hot and cold phases, inhomogeneously dis- 
tributed metals in the ICM (Fabian et al. 2001, 2002), AGN heat- 
ing by jets (Briiggen & Kaiser 2002) and sound waves (Fabian 
et al. 2003), thermal conduction (Voigt et al. 2002), and a signif- 
icant relativistic cosmic ray component frozen into the thermal gas 
(Cen 2005). 

Secondly, using more sensitive receivers, Edge (2001) has de- 
tected CO emission in 16 cooling flow central galaxies, consis- 



tent with 10 J - 10 11 -' M of H 2 at 40 K (see also Salome & 
Combes 2003). These are roughly the masses expected, given the 
revised cooling rates and likely ages. Interferometry shows fur- 
ther that the CO emission is localised within the central few arc- 
sec of the cluster (Edge & Frayer 2003; Salome & Combes 2004). 
The frequent occurrence of a much hotter molecular gas compo- 
nent in these systems has also been established, via near-infrared 
surveys of the ro-vibrational H2 lines. Building on some earlier 
work (Elston & Maloney 1992; Jaffe & Bremer 1997; Falcke et 
al. 1998; Jaffe, Bremer & van der Werf 2001), Edge et al. (2002) 
performed an H+K band spectroscopic survey of 32 line-luminous 
central cluster galaxies and showed that the H2 emission correlates 
in strength with the CO and Ha emission. Analysis by Wilman et 
al. (2002) showed that the lower-lying H2 lines are thermally ex- 
cited in dense gas (n > 10 5 cm~ 3 ) at T ~ 2000K. The inferred 
gas pressure (nT) thus exceeds 10 8 cm -3 K, some 2-3 orders of 
magnitude higher than that in either the diffuse X-ray or optical 
emission line gas. There must be a close physical connection be- 
tween the H2 and optical emission line clouds, since the Fb/Ha 
ratio is observed to be constant over 2 orders of magnitude in Ha 
luminosity. Narrow-band HST NICMOS imaging of three central 
cluster galaxies by Donahue et al. (2000) had earlier shown that the 
hot H2 and hydrogen recombination emission lines have similar 
morphologies. 

Following on from our spectroscopic survey in Edge et 
al. (2002), we have begun detailed studies of individual objects to 
analyse the relationship between the various ionized and molecu- 
lar gas systems. Here we present H and K-band area spectroscopy 
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of the centre of NGC 1275 using the UIST Integral Field Unit on 
the United Kingdom Infrared Telescope (UKIRT). From narrow- 
band HST observations of this target, Donahue et al. (2000) found 
that the bulk of the ro-vibrational emission originates in the nu- 
clear source, with some diffuse emission lying up to a kilopar- 
sec away. With the UIST IFU we can map this emission, in- 
fer its excitation state and explore its possible connection to the 
1.2 kpc radius ring of CO-emission found by Inoue et al. (1996). 
Area spectroscopy of NGC 1275 was published earlier by Krabbe 
et al. (2000) but at significantly lower spatial and spectral resolution 
and non-simultaneously in the H and K-bands. 

The structure of the results section of this paper is as follows. 
Firstly, we analyse the emission line morphology in the central re- 
gion and demonstrate that whilst the AGN continuum is (as ex- 
pected) unresolved, the H2 emission appears to be concentrated 
around ~ 50 pc from the nucleus. Thereafter, we show that con- 
sideration of the excitation and kinematic properties of the H2 are 
consistent with this interpretation, and suggest that the gas is in a 
disk-like structure spread over radii ~ 35 — 70 pc. We offer fur- 
ther support for the model by deriving a dynamical estimate of the 
mass of the central supermassive black hole, and find that it is con- 
sistent with the best non-dynamical estimates in the literature. We 
then discuss the possible connection between this H2 disk and the 
CO emission on larger scales. 

With Ho = 70 kms" 1 Mpc" 1 , the spatial scale in 
NGC 1275 (z=0.0176) is 358 pc per arcsec. 



2 OBSERVATIONS, DATA REDUCTION AND ANALYSIS 

The observations were performed on the nights of 2003 Septem- 
ber 24 and 26 using the UIST IFU on UKIRT. The IFU provides a 
6.5 x 3.4 " field of view and works on the image slicing principle, 
with the focal plane being covered by 14 'slits' which are aligned 
parallel to the long-axis of the field of view and reformatted into 
a single long slit at the entrance to the dispersion unit. The slits 
are 0.24" wide and the pixel size along the slits is 0.12". We per- 
formed observations with two separate grisms. The first observa- 
tion was taken with the short-K grism spanning 2.05-2.25 /im with 
a spectral resolution of ~ 3600, and we observed the target with 
500 second exposures for a total of 42 minutes on source (with the 
usual object-sky-sky-object nodding pattern). A second observa- 
tion was obtained with the HK grism (resolution ~ 900), spanning 
1.45-2.5/im, with 240 second exposures for a total of 52 minutes 
on source. The long-axis of the IFU was oriented in an east-west 
direction for both observations. The observations were taken under 
photometric conditions with atmospheric seeing in the range 0.3- 
0.4", which is well-matched to the spatial sampling provided by 
UIST. For convenience, the two datasets will hereafter be referred 
to as the 'short-K' and 'HK' data, respectively. 

Data reduction was performed offline using the UIST-specific 
recipes within the ORAC-DR software (Cavanagh et al. 2003) and 
comprised the usual steps of flat-fielding, sky-subtraction, wave- 
length calibration, division by an atmospheric standard and flux 
calibration. The resulting data cubes contain 14 x 54 spatial ele- 
ments for each of the 1024 spectral channels. They were cleaned 
of cosmetic defects such as hot pixels and cosmic rays by use of 
a median filtering process, whereby a two-dimensional kernel was 
passed through the cube and at each point pixels deviant from the 
median value by more than a certain threshold were replaced by 
the median. The filter was applied five times, during the course of 
which the replacement threshold was gradually lowered from 5<r 



to 2a. Subsequent cube manipulation and spectral analysis were 
performed with a combination of IDL, IRAF and QDP/PGPLOT. 



3 RESULTS 

3.1 Emission line morphologies 

We began by deriving from the HK data cube emission line maps 
of H 2 v=l-0 S(l), [FeII]A1.644 and Paa, as shown in Fig.Q They 
were formed by taking slices through the data cube at the wave- 
lengths of interest and subtracting the underlying continuum by us- 
ing line-free regions either side of the line and interpolating the 
continuum linearly across the line. The detailed spectral decompo- 
sition of the continuum by Krabbe et al. (2000) shows that there 
are no stellar absorption features in the vicinity of these emis- 
sion lines, and that our continuum subtraction procedure is there- 
fore reliable. The maps show that the emission in all three lines 
is heavily dominated by the nucleus, and although there is notice- 
ably more extended emission to the west of it than to the east, 
the 'filament' of H2 is the only feature that can be readily iso- 
lated. Its appearance in the narrow-band H2 image in Donahue 
et al. (2000) attests to its reality. The H2 v=l-0 S(l) luminosities 
of this filament and the 2 x 2" nuclear region are 1.8 x 10 39 
and 3.0 x 10 40 ergs -1 , respectively. A comparison of our mea- 
surement of the H2 v=l-0 S(l) flux in the 2 x 2" nuclear region 
(43.0 ± 0.9 x 10 -15 erg cm" 2 s" 1 ) with that measured in the cen- 
tral 3" by Krabbe et al. (2000) (41.1 ±2. 5 x 10" 15 erg cm" 2 s" 1 ), 
demonstrates that there has been no variation in the 8.7 years be- 
tween the two observations. 

The fine pixel scale of the UIST IFU, coupled with the excel- 
lent seeing during the observations, allows us to examine the inten- 
sity profiles of the emission within 1" of the nucleus and to con- 
strain the spatial locations of the emitting regions. Almost 80 and 
90 per cent, respectively, of the line and continuum emission in this 
central region falls in slits 9 and 10, the bulk of it in slit 9, which we 
refer to hereafter as the 'peak slit' . In Fig.|2|we show intensity pro- 
files along the peak slit for several emission lines, with the fluxes 
derived from single gaussian fits to the HK data. The pointspread 
function (PSF) for this slit, as derived from a reconstructed white- 
light image of a standard star, is also indicated. It can be seen that 
Paa is unresolved, whereas the H2 and [Fell] emission are both 
extended with respect to the PSF. 

A closer look at the H2 profile is shown in Fig. [5] where the 
fluxes were derived from double gaussian fits to the profiles from 
the higher spectral resolution, short-K, data cube. They show that 
the H2 profile has a FWHM = 0.6", corresponding to a PSF- 
deconvolved FWHM = 0.4". Indeed, Fig. [3] also demonstrates 
(but does not of itself prove) that the bulk of the H2 emission 
out to radii ~ 0.4" can be accounted for by a superposition of 
2 point-sources of emission, each located at a radius of approxi- 
mately 0. 15" from the continuum peak (i.e. the nucleus), which lies 
between the central two pixels. In contrast, Fig.|3|also demonstrates 
than the continuum beneath the line - dominated by AGN emission 
in this central region - is spatially unresolved. Comparison of the 
continuum fluxes in slits 9 and 10 (not shown) implies that the see- 
ing disk of the unresolved point source is offset by ~ 0.06" from 
the centre of slit 9, towards slit 10. This reflects the precision with 
which a source can be acquired within the IFU aperture and hence 
also limits our determination of the PSF from a single standard star 
observation. 

To summarise our analysis of the emission line morphologies, 
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Figure 1. Contour plots of the H2 v=l-0 S(l), [FeII]A1.644 and Paa emission derived from continuum-subtracted narrow-band cuts on the 
data cube obtained with the HK grism. The narrow-band images were smoothed with 0.48" square top-hat filters prior to contouring. The 
lowest contour levels are set at 2a of the background in the lower part of each image, and have values of 3.1 x 10 -20 , 5.2 x 10 -20 and 
1.1 x 10 -19 W m -2 , for H2 v=l-0 S(l), [FeII]A1.644 and Paa, respectively. Successive contour levels increase by \/2. The boxes used for 
extraction of spectra for the nucleus and a filamentary structure are indicated on the H2 v=l-0 S(l) map. The orientation on all three panels 
is as indicated. 
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Figure 2. Variation in emission line flux along the IFU slit passing 
through the nucleus, for Paa (crosses), [FeII]A1.644 (diamonds), 
H2 v=l-0 S(l) (triangles) and S(2) (asterisks). The lines show 
scaled PSFs for this slit, determined from a standard star obser- 
vation. Line fluxes were derived from gaussian fits to the profiles in 
the HK data cube. 

we adopt as a working hypothesis a model in which the bulk of the 
H2 emission arises in an extended structure located ~ 0.15" (50 pc) 
either side of the nucleus along an east- west axis. Further support 
for this model comes from an analysis of both the excitation and 
kinematics of the H2 emission, which we pursue below. 

3.2 H 2 excitation 

In Fig.|4|we show HK spectra of the nucleus and the filament re- 
ferred to in Fig.Q Also shown is the spectrum of the extended 
emission within the 2 x 2" nuclear box, formed by subtracting a 



scaled version of the spectrum of the central 9 pixels (scaled using a 
PSF observation of a standard star to take account of light from the 
point source scattered into other pixels). Although it falls in a re- 
gion of bad atmospheric cancellation, Paa appears to be very weak 
in the circumnuclear spectrum, confirming the result of Fig.|2|that 
it is spatially unresolved. 

To examine the excitation mechanism for the H2, we begin 
by comparing the observed H2 line ratios with the predictions for 
thermal (collisional) excitation, which dominates in molecular gas 
with density tit > 10 5 cm -3 (the critical density of these transi- 
tions), heated to a few thousand Kelvin by UV or X-radiation, or 
in shocks. Under these conditions, the occupation numbers of the 
excited ro- vibrational levels of the H2 molecule will be in thermal 
equilibrium at a temperature T cx equal to the kinetic temperature of 
the gas. Therefore, for a given object, the flux f] in the ith emission 
line will satisfy the relation Zog(FiAi/Ai5fi) = constant — Ti/T ex , 
where Ai is the wavelength of the line, and Ai, TS and gi are the 
spontaneous emission coefficient, energy (expressed as a tempera- 
ture) and statistical weight of the upper level of the transition (the 
latter assumes an ortho:para H2 abundance ratio of 3: 1, appropriate 
where collisions dominate). 

Such plots were constructed for the nucleus, the circumnuclear 
region and the filament, and are shown in Fig. [5] For the nucleus, 
we see that the low-lying transitions are completely thermalised, 
with an excitation temperature of 1360 ± 50 K, implying a gas 
density above 10 5 cm -3 . The excess flux in the higher-excitation 
transitions may have one of three possible origins: (a) it could be 
from a much hotter thermalised gas with T ox = 3100 ± 500 K (as 
deduced from a fit to the high-excitation lines only), close to the 
H2 dissociation temperature of ~ 4000 K (fitting instead a two- 
temperature model with the excitation temperature of the hotter 
component fixed at 3100 K implies that the cooler component has a 
temperature of 800 ± 50 K) (b) it could be non-thermal fluorescent 
emission in low density gas; (c) it could be non-thermal emission 
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due to excitation of the molecule by secondary electrons deep in the 
cloud. We confirm the finding of Krabbe et al. (2000) that on the ba- 
sis of model (a) the v=2-l S(3) transition with 71 = 13890 K is un- 
derluminous by a factor of ~ 3, explained by the fact that the upper 
level of this transition is de-populated by an accidental resonance 
with photons around the wavelength of HI Lya at 1216A (Black & 
van Dishoeck 1987). Our excitation temperature for the low excita- 
tion lines is slightly lower than that found by Krabbe et al. (2000), 
who measured 1480 ± 250 K, but our determination is more pre- 
cise due to our inclusion of several 1-0 Q series lines above 2.4/im. 
For the higher excitation lines, Krabbe et al. determined tempera- 
tures between 2600 and 2900 K. For the circumnuclear and filament 
and circumnuclear spectra, we deduce excitation temperatures of 
2200 K, as expected for X-ray or shock-heating. 

Given an excitation temperature, the v=l-0 S(l) fluxes can be 
converted to masses of hot molecular hydrogen in the following 
manner. For optically thin emission (which is a valid assumption 
for these transitions), the H2 column density, N(H2), is related to 
the observed intensity of the v=l-0 S(l) emission, I, by the ex- 
pression N(H 2 ) = 4%1/fAhv, where A = 3.47 x lO^s" 1 , 
hv — 9. 37x 10 -13 erg, and / is the fraction of H2 molecules in the 
v = 1, J = 3 state leading to v=l-0 S(l) emission. Integration over 
the solid angle on the sky covered by the source implies that the 
mass of hot H2 is given by M(i?2) = 4:nL 1 _ s(i)mH 2 / 'j ' Ahv, 
where I/i_os(i) is the v=l-0 S(l) luminosity and m,H 2 is the mass 
of the H2 molecule. Computation of the ro- vibrational H2 partition 
function at various temperatures enables us to deduce the appropri- 
ate value of /, which varies as follows: 0.00031 (800 K), 0.0049 
(1360 K), 0.016 (2200 K), 0.023 (3100 K). Hence for the nucleus 
we deduce the following H 2 masses: (i) 3.9 x 10 5 M©, if all the 
v=l-0 S(l) is thermal emission at 1360 K; (ii) 1.0 x 10 6 M Q at 
800 K and 4.5 x 10 4 M Q at 3100 K, in the aforementioned two- 
temperature model. For the circumnuclear and filament spectra, the 
masses of hot H 2 at 2200 K are 3.3 x 10 4 M Q and 7.3 x 10 3 M , 
respectively. 

With reference to the work of Maloney, Hollenbach & Tie- 
lens (1996) we now assess the energetic requirements for X-ray 
heating of the molecular gas. As discussed therein, the dominant 
parameter controlling the physical conditions is -Hx/n, the ratio 
of the X-ray energy deposition rate per particle, Hx, to the gas 
density, n. The former is the integral over energy of F(E)a pe , 
where F(E) is the local photon energy flux per unit energy in- 
terval and tjpo the photoelectric cross-section per hydrogen atom. 
For an X-ray source with a power-law spectrum of photon index 
r = 2 and 1-100 keV luminosity L x = 10 44 L 44 ergs' 1 , located 
50rso pc from the gas cloud: 

/Lx ~ 3 x lQ- 21 L iiT ^N^ ergs" 1 (1) 

where Nu = 10 22 A^22 cm -2 is the equivalent neutral hydro- 
gen column density attenuating the X-ray flux; a column of at least 
10 21 cm -2 is required in order to exclude the normal UV photon- 
dominated region at the cloud surface. Applied to NGC 1275, we 
find that: 

logHx/n ~ -25.3 + log[L x /L E instcin ] - log[r 50 n 5 N 2 2]- (2) 

Here Lx is the X-ray luminosity of the source as seen by the 
cloud and I/Einstoin = 1.6 x 10 44 ergs -1 is the X-ray luminos- 
ity of the source extrapolated (in energy) from the measurement by 
the Einstein satellite in 1979 (Branduardi-Raymont et al. 1981), 
which also suggests an obscuring column N22 ~ 1; n$ is the 
gas density in units of 10 5 cm -3 , which as mentioned above must 
be unity or higher for thermal excitation. The modelling of Mal- 
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Figure 3. The diamonds show the H2 v=l-0 S(l) flux in pixels along 
the peak IFU row measured relative to the continuum nucleus (neg- 
ative offsets are to the east). The dashed lines show the contribu- 
tions from two point-sources of emission at ±0.12", and the solid 
line their sum. The latter accounts for the bulk of the emission in 
the central 0.6". Line fluxes were derived from double gaussian fits 
to the profiles in the short-K data (see section 3.3). The asterisks 
connected with the dotted line show (with arbitrary normalization) 
the continuum level below the H2 v=l-0 S(l) line: it is clearly less 
extended than the line emission and spatially unresolved. 

oney et al. (1996) demonstrates that the bulk of the thermally ex- 
cited H2 emission is produced over a fairly narrow range around 
logHx/n — —26, with the H2 emissivity dropping by at least a 
factor of 3 when logHx/n is more than ±0.3 dex from the emis- 
sivity peak. Equation (2) therefore demonstrates that X-ray heating 
is an energetically feasible excitation mechanism for the observed 
thermal H2 emission. It is, however, difficult to make a more pre- 
cise statement because of the uncertainy over the appropriate value 
of Lx, since the light travel time from the nucleus to the clouds is 
~ 160 yr, but observations of the source over just the last 30 years 
show that Lx has been far from constant. From the early- 1970s 
to the mid-1990s the X-ray luminosity of the central source has 
dropped steadily by a factor of 20 (see e.g. Levinson et al. 1995 for 
a summary of the observational evidence). It is, however, clear that 
unless the AGN were much weaker in the past or the gas signifi- 
cantly denser than ns ~ 1, the hot H2 could not exist much within 
the observed 50 pc radius. Similarly, equation (2) also dictates that, 
all other parameters being equal, thermally excited H2 emission is 
produced in fairly narrow range of radii, within ±0.15 dex of the 
peak radius, ~ 35 — 70 pc in this case. 

In concluding this subsection, we point out that the results 
of Maloney et al. (1996) for X-ray heating assume an equi- 
librium state, which is established on a time-scale of i cq ~ 
750L7 4 1 r 2 oyr (set by the ionization rate of hydrogen, which has 
the slowest chemistry to equilibriate). Therefore, if the central X- 
ray source is varying significantly on timescales of t aq or less - and 
observations over the last 3 decades suggest that it is - considera- 
tion of non-equilibrium chemistry may be necessary. 



3.3 H2 kinematics: evidence for a rotationally-supported 
structure 

We turn our attention now to the kinematics of the molecular gas 
in the vicinity of the nucleus, as derived from the short-K data. The 
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Figure 4. HK spectra of the nucleus, the filament and circumnuclear region in NGC 1 275 . The boxes used for the extraction of the nucleus and filament spectra 
are shown in Fig.Q the spectrum of the circumnuclear region was obtained by subtracting from the 2 X 2" nuclear box an appropriately scaled version of the 
spectrum of the central 9 pixels, such that only the extended emission in this region remains. Note the broad base to the Paa line on nucleus, consistent with 
the Seyfert 1 optical classification of NGC 1275. 



most striking feature, visible in the raw data, is the strong veloc- 
ity shear in the peak position of the H2 v=l-0 S(l) line moving 
across the nucleus. To demonstrate this, we show in Fig. [6] line 
profiles on a pixel-by-pixel basis along the peak slit of the IFU, 
passing through the nucleus. The kinematics were quantified by fit- 
ting double gaussian profiles to the lines, with the widths of the two 
components constrained to be the same. The fitting of such profiles 
does not imply that gas of both velocities is present at each physical 
location in the galaxy; rather, the weaker of the two components is 
simply scattered light from gas on the other side of the nucleus, in 
line with the model of Fig. [5] in which the bulk of the H2 is con- 
centrated in two approximately point-like components 0.15" either 
side of the nucleus. Further support for the latter model is provided 
by Fig. Q which shows the fluxes of the fitted red and blue com- 
ponents in the peak slit (slit 9) and also in slit 10. We also show 
rotation curves for the central four slits (i.e. slits 8-11); for slits 9 
and 10, the velocity being that of the brighter of the two gaussians 
at each location. For slits 8 and 11, the lines are fainter, so only a 
single gaussian can be fitted to these lines. The appearance of these 
rotation curves, with the sharp decrease in velocity amplitude in the 
slits either side of the nucleus suggests that the molecular gas may 
be in a disk-like structure, rotating about an axis oriented approxi- 
mately north-south, coincident with the radio-jets. 

What is the nature of this rotating H2 structure? It is unlikely 
to be part of an optically-thick torus of the type commonly postu- 
lated in the centres of Seyfert galaxies (e.g. Antonucci et al. 1993), 
although there is evidence for such a parsec-scale torus in NGC 



1275. Firstly, the presence of a hot dust continuum in the K-band 
implies that some dust is heated to its sublimation temperature and 
thus lies as close as 0.05 pc to the nucleus (Krabbe et al. 2000). 
Secondly, gas in a parsec-scale torus perpendicular to the radio jet 
axis is a plausible candidate for producing the observed free-free 
absorption towards the milliarcsec-scale northern radio counter-jet 
in NGC 1275 (Levinson, Laor & Vermeulen 1995); with reference 
to the latter's models, constraints since obtained on the black hole 
mass in NGC 1275 (see subsection 3.5) suggest that the torus is 
comprised of individual clumps of dense gas (n 2> 10 4 cm -3 ) 
and is matter-bounded in the radial direction. At larger radii, NGC 
1275 has a complex dust distribution extending out to 17 kpc, in 
contrast to the majority of FRI radio sources in which the dust is 
settled in disks with radii < 2.5 kpc (de Koff et al. 2000) 

A more likely possibility is that the hot H2 is part of a disk, 
from which we only see emission from a ring covering a narrow 
range of radii around 50 pc, due to the nature of the X-ray heat- 
ing by the central source (see section 3.2). Would such a disk be 
stable? We consider first a gas disk, in vertical hydrostatic equi- 
librium in the gravitational field of a central black hole of mass 
M. To ensure that the disk is stable against self-gravity we re- 
quire E < |pf, where E is its surface mass density and h its 
exponential scale-height (this is essentially the Toomre (1964) sta- 
bility criterion). There are two approaches to evaluating this cri- 
terion, leading ultimately to the same conclusion: (i) integrating 
vertically over the disk density profile, n = noexp(—z 2 /h 2 ), 
we find E = 2^/TTnohm p . Hence the stability criterion becomes 
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r 3 < - — — — t=, The mid-plane density, nn must be at least 

10 5 cm -3 for thermal H2 excitation, and taking M = 3.4 x 
10 8 M as the black hole mass (see section 3.5), this translates 
into r < 21 pc i.e. less than the observed radius of the H2 distri- 
bution, (ii) Assuming that the mass of hot nuclear H2 deduced in 
section 3.2 (3.9 x 10 Mq) is spread over an annulus from 35- 
70pc we deduce E = 0.07 Kg m~ 2 , so re-casting the stability 
criterion with r — 50 pc we require h/r > 10~ 3 . However, from 
above E = 2y / 7rno/im p , so requiring no > 10 5 cm -3 we deduce 
that h/r < 8 x 10 , once again inconsistent with the stability cri- 
terion. We therefore conclude that a smooth, stable, gaseous disk of 
the observed radius and mass density cannot exist. The propensity 
for lOOpc-scale nuclear gas disks to be subject to gravitational in- 
stabilities was highlighted recently by Tan & Blackman (2004) in 
the context of accretion onto the nuclei of giant elliptical galaxies. 
The disk can, however, be stabilised if the gas is located in dense 
clouds, as opposed to being smoothly distributed. In this case, the 
Toomre stability criterion can be cast as ^ > 1.68, where a is 
the 1-d velocity dispersion of the clouds and Q = y/ GM/r s the 
orbital angular velocity. Assuming E = 0.07 Kg m~ 2 as before, 
we hence require a > 70 ms -1 for stability, which is not at all 
restrictive. 

From the short-K data we also investigated the kinematics of 
the H2 filament. With respect to the assumed velocity zero-point in 



Fig.|7| the H2 velocity of this gas is —10 ± 14kms -1 , with a de- 
convolved linewidth of 190 km s -1 FWHM. Hence, this material 
is unlikely to be part of the putative disk. The wavelength range of 
the short-K data does not contain [Fell] or Paa so we are unable to 
study the kinematics of these lines in comparable detail. 

In this subsection we have demonstrated that the H2 is plausi- 
bly in a disk-like structure, composed of individual clumps of dense 
gas since a smooth gas disk of the observed dimensions would be 
gravitationally unstable. Because the points in the rotation curves 
are not mutually independent - the pixel separation along each slit 
is 0.12" whilst the seeing is ~ 0.4" - we have not attempted to fit 
a full two-dimensional disk model to the data, since this would re- 
quire a much better knowledge of the PSF than we have. Adaptive 
optics observations are needed to conclusively test the disk model 
in two dimensions. Nevertheless, in the next section we assume that 
the disk model is valid and make an estimate of the black hole mass. 



3.4 Black hole mass estimation 

The computation of dynamical black hole mass estimates for galac- 
tic nuclei has reached a considerable level of complexity, requiring 
in general observations of the full two dimensional velocity and 
velocity dispersion fields, knowledge of the distributed stellar mass 
and dynamical modelling. The resulting model of the gas kinemat- 
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ics must then be projected on to the plane of the sky, convolved with 
the instrumental response and fitted to the observational data. Here 
we circumvent all of this complexity and compute a simple estimate 
of the mass enclosed within the observed radius of the molecular 
gas. Thereafter, we use an archival HST NICMOS image of NGC 
1275 to estimate the stellar mass contribution in this region. The 
resulting black hole mass is then found to compare favourably with 
other non-dynamical estimates. 

From the rotation curve for the nuclear slit (slit 9) shown in 
Fig.Q we see that the H2 v=l-0 S(l) velocity jumps by approxi- 
mately 240 kms" 1 across the nucleus. If the perpendicular to the 
disk is inclined at an angle i to the plane of the sky the velocity 
difference increases to AV = 240/sin i kms -1 . Our best esti- 
mate of i derives from VLBI radio observations of the parsec-scale 
jets, which suggest that the radio jets are inclined to the line of 
sight at an angle of 30-55 degrees, with the southern jet approach- 
ing us (Walker, Romney & Benson 1994). On the assumption that 
the molecular gas rotation axis coincides with the radio jet axis we 
take i = 45 degrees. Within the inferred molecular gas radius of 
r — 50 pc, the dynamical mass is given by: 



M (total) 



G 



(3) 



with v = AV/2 = 120/sin i kms x . The result is 
M (total) = 3.4 x 10 8 Mq. Taking account of the aforementioned 
uncertainty in the inclination of the system introduces a systematic 
uncertainty of ±0.18 dex on this mass. 

For the next stage of the calculation we used an image of 
NGC 1275 to infer the stellar mass distribution. The data in ques- 
tion are a 640 second exposure HST NICMOS image taken in the 
F160W filter (roughly corresponding to the H-band) on 03/16/98 
(dataset number N3ZB1R010). The standard pipeline-processed 
dataset was obtained from the archive. Using the IRAF STSDAS 
task ellipse, elliptical isophotes were fitted and the resulting sur- 
face brightness profile is shown in Fig. [8] Since the eccentricity of 
the isophotes is typically small (e ~ 0.2) we chose to fit the profile 
with a spherically symmetric modified Hubble Law: 



I(R) 



Io 



1 + (R/Ro) 



(4) 



which deprojects analytically to the following spatial luminos- 
ity density: 



j(r) = 



Jo 



[1 + (r/7? ) 2 ] 3 / 2 ' 



(5) 



Fitting to the data at radii > 0.5" (beyond the region 
dominated by the AGN power-law and re-radiated dust emis- 
sion; see Krabbe et al. 2000 for a spectral decomposition of 
the continuum) we deduce Ro = 2.3" and Io = 1.12 x 
10 erg cm s A arcsec . Although there are some sys- 
tematic deviations from the model at R > 3" it is certainly good 
to within a factor of 2. When deprojected, the implied central stel- 
lar mass density is 6.3(M/L)_h Mq pc -3 , where (M/L)h is the 
mass-to-light ratio of the stellar population in solar units. Values 
of (M/L)h ~ 1 are expected for stellar populations with age 
~ 10 Gyr for a wide range of power-law and exponential stel- 
lar initial mass functions and metallicities (Salasnich et al. 2000); 
also, for a given population, (M/L)h will vary in the approxi- 
mate range 0.05-3 as the population evolves, and will in general 
increase with age. Assuming, therefore, that (M/L)h = 1, the 
implied stellar mass within the 50 pc radius of the molecular gas 
distribution is 3.7 x 10 6 Mq, which is ~ 1 per cent of the above 



figure for Al (total). Hence our estimate of the black hole mass is 
M(BH) = 3.4 x 10* M Q , ±0.18 dex. 

It is interesting to compare this figure with other estimates of 
the black hole mass in NGC 1275. Bettoni et al. (2003) used the 
observed central stellar velocity dispersion (o — 250 km s -1 ) in 
conjunction with the established M(BH) — a correlations (Geb- 
hardt et al. 2000; Ferrarese & Merritt 2000) to deduce M(BH) = 
4.1 x 10 s M Q , ±0.4 dex. Using the M(BH) - L Bu i g c relation 
(Magorrian et al. 1998), in which the scatter is arguably some- 
what higher, they calculate M(BH) = 1.4 x 10 9 Mq, ±0.4 dex. 
For comparison, dynamical estimates of the black hole mass are 
available for two other nearby central cluster galaxies, Cygnus A 
(M(BH) = 2.5 ± 0.7 x 10 9 M G ; Tadhunter et al. 2003) and M87 
(M(BH) = 3.2±0.9x 10 9 M ; Marconi et al. 1997). NGC 1275 
and M87 are both FRI radio sources with virtually the same mod- 
est radio power (logP(178MHz) = 24.6 W nT 2 sr _1 ), but the 
large difference in their black hole masses is consistent with the 
scatter at low power in the large sample of z ~ 0.5 radio galax- 
ies in McLure et al. (2004). Cygnus A, in contrast, is a powerful 
FRII radio source (logP(l78MHz) = 27.7 W m" 2 sr _1 ) with a 
high-excitation nuclear spectrum, and falls at the high end of the 
M(B J/)-radio power correlation for such objects in the McLure 
et al. sample. We conclude that our dynamical estimate of the black 
hole mass in NGC 1275 is consistent with the most reliable existing 
estimate and that this value is itself in line with expectations based 
on its radio properties. 



3.5 Relationship between the hot H2 and the large-scale 
distribution of CO emission 

Interferometry of the CO J=l-0 line by Inoue et al. (1996) sug- 
gests that some 3 x 10 10 Mq of cool H2 exists within the central 
arcmin of NGC 1275, in the form of a plume extending from the 
nucleus to lOkpc west of it. Within the vicinity of the nucleus, 
6 x 10 9 Mq of this material is confined within a ring-like structure 
which appears on the sky as two peaks: peak 1, 1.2 kpc west of the 
nucleus, and peak 2 the same distance to the south east. The gas ap- 
pears to be rotating about the nucleus with line-of-sight velocities 
±150 kms -1 , in the same sense as the H2 we have discovered on 
much smaller scales, suggesting a possible dynamical connection 
between the two gas systems. Inoue et al. (1996) suggested that the 
CO ring could be due to gas trapped at the inner Lindblad resonance 
(ILR) of the galaxy potential; from poorer spatial and spectral res- 
olution near-infrared observations they suggested that the central 
peak of hot H2 could be due to shock emission in a turbulent distri- 
bution of molecular clumps funnelled inwards from the CO ring. In 
contrast, we find that the small-scale hot H2 is in a rotationally sup- 
ported structure and that the nuclear radiation field is the principle 
source of excitation. 

To further investigate the connection between the CO ring and 
the hot H2 disk, it would be of interest to search for H2 emission co- 
incident with the former (which falls outside out UIST IFU field of 
view). There is no H2 emission at this position in the narrow-band 
HST NICMOS H 2 image of Donahue et al. (2000), suggesting that 
any associated emission is fainter than the H2 'filament' (which 
they do detect; see our Fig.0. This filament may represent (pos- 
sibly shocked) material being transported from the CO ring to the 
H2 disk: it lies at the systematic velocity (see section 3.3) and has a 
thermal H2 excitation characteristic of shock emission (see section 
3.2). 

We also note that the absolute velocity of CO "peak 1" 
1.2 kpc west of the nucleus is virtually identical to that of the 



8 R.J. Wilman et al. 



100.00 



o 

0) 




CO 

CO 



0.01 I 

0.1 1.0 10.0 

Radius(arcsec) 

Figure 8. The surface brightness of the inner regions of NGC 1275 
as function of the radius of the semi-major axis of the ellipti- 
cal isophotes, derived from an HST NICMOS observation in the 
F160W filter (corresponding closely to the H-band). The line shows 
the best fit modified Hubble profile at radii > 0.5", with the upturn 
in the data below this radius being caused by non-stellar emission 
associated with the AGN. 



have large masses of CO localised in the central few tens of kpc and 
smaller masses of hot H2 on even more compact scales. Regardless 
of the precise origin of the cold gas (be it from the cooling flow 
itself, or an interacting galaxy), the CO-H2 connection found here 
for NGC 1275 may be universally important for fuelling AGN in 
the central galaxies of cooling flows, and hence for the regulation 
of heating and cooling in such cluster cores. 

These results demonstrate the utility of the UIST IFU for high 
spatial and spectral resolution work of this kind, including black 
hole mass measurements, and promise much for the advent of sim- 
ilar instruments on 8-metre telescopes with adaptive optics. When 
coupled with sensitive CO interferometry (e.g. from ALMA) it 
should be possible to extend these detailed studies to more distant 
clusters. 
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hot H2 at a radius of 50 pc on the same side of the nucleus: on 
the slit 9 rotation curve of Fig. CO "peak 1" would lie at an 
offset of +3.5" with velocity — 125 kms -1 . In the most naive 
of dynamical models this would imply that M(< r)/r is the 
same at radii of 50 pc and 1.2 kpc. Integrating equation (5), the 
enclosed stellar masses within these radii are 3.7 x 10 6 (MIL) 11 
and 1.5 x 10 10 (M/L)h, respectively. Equality in M(< r)/r can 
be achieved by adding a central black hole of mass M(BH) = 
6.5 x 10 8 (M/L) H M q . This shows that the black hole mass de- 
duced from consideration of the hot H2 dynamics alone is not in- 
consistent with the assumption that (M/L)h is close to unity. 



4 CONCLUSIONS 

For the first time, we have spatially and kinematically resolved the 
hot H2 in the circumnuclear region of NGC 1275 and found that 
it could be concentrated in a rotationally-supported structure in a 
narrow range of radius around 50 pc from the nucleus. A smaller 
amount of extended H2 emission exists within a few hundred pc of 
the nucleus. X-ray heating by the active nucleus has been demon- 
strated to be a viable heating mechanism for the bulk of the ther- 
mally excited H2 emission, and also explains why the emission is 
dominated by gas within a very narrow range of radii. Thanks to 
the fine pixel scale of the UIST IFU and excellent seeing, we were 
able to observe a sharp shift in the velocity of the H2 v=l-0 S(l) 
across the nucleus and to make a simple dynamical measurement 
of the black hole mass, of 3.4 x 10 s M©. The principal uncertainty 
in the mass arises from the uncertainty in the true inclination of 
the molecular gas disk to the plane of the sky. Nevertheless, the 
measured value is within 20 per cent of that inferred by Bettoni et 
al. (2003) from the M(BH) - a relation. 

The H2 structure is likely to be an extension to smaller spatial 
scales of the coaxial 1.2kpc-radius ring of CO emission found by 
Inoue et al. (1996), which is itself the terminus of a lOkpc-long 
plume of CO. This suggests that we may be directly witnessing the 
delivery of fuel from the galactic scale to the AGN itself. As noted 
in section 1, increasing numbers of cooling flow clusters appear to 
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Figure 6. Emission line profiles of the H2 v=l-0 S(l) line along the peak IFU slit in the short-K data, which passes east-west through the nucleus, starting 
from pixel (9,27) in the east to (9,32) in the west. Adjacent pixels are separated by 0.12" and the nucleus itself is assumed to lie mid-way between pixels 
(9,29) and (9,30) where the continuum peaks. Fits to the line with double gaussian profiles are also shown, and note the sharp shift in the velocity of the line 
peak across the nucleus. 
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Figure 7. Upper left: the diamonds show the fluxes of the red and blue gaussian components fitted to the H2 V=l-0 S(l) lines in slit 9 (see Fig.|6|, and the thick 
solid line their sum. The dashed lines are PSFs. Upper right: as at upper left but for slit 10 (in both these slits, the blue component is the one which dominates 
at positive offsets) . Centre left: the velocity of the brighter of the two H2 v=l-0 S(l) gaussian line components as a function of position for slit 9. Centre right: 
as at centre left but for slit 10. Lower panels: rotation curves for slits 8 and 1 1, as derived from single gaussian fits to H2 V=l-0 S(l). The abscissa on all plots 
is the offset along the slit relative to pixel position 29.5 (see Fig.|S|, with negative offsets corresponding to points east of the nucleus. All velocity error bars 
are lcr. 



